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This paper discusses the generation of a high precision DEM (Digital Elevation Model) based on high 
density airborne LiDAR (Light Detection and Ranging) data for an interdisciplinary landscape archaeolog¬ 
ical study concerning the settlement history and environment in Sandy Flanders, a region to the north of 
Ghent (Belgium). The objective was to create a detailed topographical surface free of artificial features 
and topographical artefacts, in the form of a DEM, visualizing the natural and current topography through 
the implementation of true ground points only. The semi-automatical removal of these features and arte¬ 
facts was based on topographical vector data, visual interpretations and slope analysis. Ultimately two 
DEM’s were constructed (1) a TIN (Triangulated Irregular Network) model, whereby the inherent large 
file format restricts the usability to large scale and (2) a grid model which can be used for small-, 
medium- and large-scale applications. Both datasets were used as an image that is interpreted using 
ancillary data from historical sources. Its usefulness is illustrated in a case of field pattern and microfield 
topography. Starting from this DEM, the approach of this landscape historical study is mainly retrogres¬ 
sive, i.e. starting from the landscape structures and elements that are still present in the contemporary 
landscape and moving into the past. 

© 2010 Published by Elsevier Ltd. 


1. Introduction 

Airborne LiDAR data for topographic data collection is available 
since the 1980s (Ackermann, 1999; Drosos & Farmakis, 2006; Liu, 

2008) . This technology is frequently used in a broad range of appli¬ 
cations, such as geomorphological research (Eeckhaut, Verstraeten, 
& Poesen, 2007; Kasai, Ikeda, Asahina, & Fujisawa, 2009), coastal 
zone monitoring (Fernandes, Utkin, Lavrov, & Vilar, 2004; Lohr, 
1998; Mutlu, Popescu, & Zhao, 2008; Saye, Wal, Pye, & Blott, 
2005; Utkin, Lavrov, Costa, Simoes, & Vilar, 2002; Zhou & Xie, 

2009) , forest fire management (Fernandes et al., 2004; Mutlu et 
al., 2008; Utkin et al., 2002) and infrastructural and environmental 
projects (Challis, Kokalj, Kincey, Moscrop, & Howard, 2008; Wehr & 
Lohr, 1999a). More significant applications for this technology are 

Abbreviations: DEM, Digital Elevation Model; LiDAR, Light Detection and 
Ranging; TIN, Triangulated Irregular Network. 
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found in archaeological surveys (Bewley, Crutchley, & Shell, 2005; 
Carey, Brown, Challis, Howard, & Cooper, 2006; Challis, 2006; 
Devereux, Amable, & Crow, 2008; EnglishJTeritage; Gallagher & 
Josephs, 2008; Harmon, Leone, Prince, & Snyder, 2006; Powlesland 
et al., 2006) and to support archaeological management programs 
(Bewley et al., 2005). 

LiDAR is an active remote sensing technique (Mutlu et al., 2008) 
which provides sub-randomly distributed three-dimensional ter¬ 
rain point data (Axelsson, 1999; Drosos & Farmakis, 2006; Liu, 
2008; Lloyd & Atkinson, 2005). It uses a laser pulse to measure 
the terrain elevation derived from the time between emission 
and reception of reflected pulses (Drosos & Farmakis, 2006; Liu, 
2008; Lloyd & Atkinson, 2005). Through laser altimetry, the x-, 
y- and z-coordinates are registered. The x- and y-coordinates rep¬ 
resent the planimetric position and the z-coordinate the altimetric 
position (Drosos & Farmakis, 2006). The final data is generally pro¬ 
vided in a comma-delimited ASCII format (Drosos & Farmakis, 
2006; Liu, 2008). Because only coordinates are provided and no 
object information is given, the interpretability of the raw data is 
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limited (Axelsson, 1999). Moreover, LiDAR data acquisition is 
expensive and consists of enormous volumes of data (Axelsson, 
1999; Challis et al., 2008; Liu, 2008). However, the exceptional pla- 
nimetric accuracy of centimeter level, which is the main advan¬ 
tage, and the regular scan pattern with high density allow the 
production of a high quality Digital Elevation Model (DEM) reduc¬ 
ing data volume significantly and making the elevation data more 
workable (Axelsson, 1999; Drosos & Farmakis, 2006; Liu, 2008; 
Lohr, 1998). Because of this, airborne LiDAR is replacing gradually 
conventional stereophotogrammetric methods (Drosos & Farmakis, 
2006; Liu, 2008). An additional advantage of LiDAR altimetry com¬ 
pared to photogrammetry is that it also penetrates vegetation, by 
the multiple echo characteristics and allows acquiring ground level 
points trough vegetation cover (Ackermann, 1999; Axelsson, 1999; 
Drosos & Farmakis, 2006; Liu, 2008; Wehr & Lohr, 1999b; Zhou, 
Song, Simmers, & Cheng, 2004). In European type coniferous and 
deciduous forests, penetration rates of 20 till 40% can be expected 
and up to nearly 70% in deciduous forests in winter time (Acker¬ 
mann, 1999). The fast and slightly higher automated LiDAR data 
processing is moreover much cheaper compared to vertical aerial 
photography (Ackermann, 1999; Axelsson, 1999; Challis et al., 
2008; Liu, 2008). For details concerning the technical principles 
of LiDAR, readers are referred to Drosos and Farmakis (2006), Wehr 
and Lohr (1999a) and Lohr (1998). 

The study area fits within the interdisciplinary landscape 
archaeological project concerning settlement history and environ¬ 
mental change in Sandy Flanders north of Ghent (Belgium). The 
overall goal of this project is to understand the spatial distribution 
of archaeological artefacts of different pre- and protohistoric peri¬ 
ods in relation to landscape and environmental changes. Princi¬ 
pally, the project studies how the palaeo-landscape determined 
the occupation and land exploitation strategies, starting from the 
Late Glacial up until the beginning of the Roman period. This re¬ 
quires an interdisciplinary, landscape-orientated approach. Several 
research processes are carried out simultaneously by different re¬ 
search groups of the Ghent University, in particular Archaeology, 
Geography, Palaeontology, Soil Management, Geology and Soil 
Science. In terms of historical landscape conditions and landscape 
change over time, this knowledge is particularly valuable for 
successful landscape protection and planning (Bender, 2009). Not 
only the historical development of the present landscape but also 
a reconstruction of past landscape conditions can be useful for 
applied historical geography (Bender, 2009). In an international 
context, the project fits into the broader framework of the Euro¬ 
pean treaty of La Valetta of the Council of Europe (1992) concern¬ 
ing archaeological heritage. The project will contribute to the 
development of scientific methods required for assessing the 
archaeological potential of areas under threat by present-day 
activities. 

The specific objectives of this paper are twofold. Firstly, trans¬ 
forming airborne LiDAR data into a high precision DEM of the re¬ 
cent and natural surface free of artefacts for landscape historical 
research and secondly the use of the DEM for the study of field pat¬ 
tern and microfield topography. 


2. Materials and methods 

2.1. Study area 

The study area covers 1400 km 2 and ranges from the city of 
Bruges in the west to the town of Temse in the east and stretches 
north of the city of Ghent to the northern state border with The 
Netherlands (Fig. 1). The area corresponds to the northern part of 
the Pleistocene valley of the Scheldt River, filled up with sandy, 
Pleistocene deposits which were remodeled during the Early 


Holocene (Heyse, 1983). Different landscape types occur within 
this generally flat area, characterized by a subtle micromorphol¬ 
ogy. This variation was an important factor in the settlement 
history as illustrated by the discontinuous spread, both chronolog¬ 
ically and geographically, of archaeological sites (Fig. 1). Important 
landscape types are a late glacial fossil lake, filled with peat and 
lake marl, cuestas of Tertiary clay outcrops, river and brook valleys, 
a massive dune complex and some smaller sandy ridges alternated 
with small depressions. Soils and environmental conditions vary 
with geomorphology as well. 

The integrative geoarchaeological project uses three scales of 
analysis. The comprehensive area (zone 1 in Fig. 1) covers the 
northern part of the Pleistocene valley. This area is not only impor¬ 
tant as a context for more detailed scientific analysis, it is also valu¬ 
able for heritage protection and management. According to the 
Flemish landscape atlas, many relics of traditional landscapes still 
exist here (Antrop & Eetvelde, 2007). The meso-level (zone 2 in 
Fig. 1) covers the hydrographical basin to which the smaller 
depression of zone 3 belongs. On this mesoscale spatial modelling 
of the development of the hydrographical network and a land eval¬ 
uation related to the archaeological sites is performed. This zone 
contains heritage landscapes that are subject to spatial planning 
procedures according to the implementation of the European Land¬ 
scape Convention (Antrop & Eetvelde, 2007). The microscale corre¬ 
sponds to a former late glacial lake and lake environment where 
peat was formed, which was an important natural resource during 
history (zone 3 in Fig. 1). The morphogenesis of this marshy 
depression is not well known yet. On this micro-level all research 
teams work together using different approaches and the most de¬ 
tailed data sources and intensive field work. 

2.2. LiDAR data and DEM generation 

The high precision airborne LiDAR measurements are distrib¬ 
uted by the Flemish GIS Agency (AGIV) and were obtained for 
the entire macro-scale study area (zone 1 in Fig. 1). The data is reg¬ 
istered between 2001 and 2004 in the scope of integral water man¬ 
agement. They consist of flight strips with swath widths varying 
between 350 and 750 m and a sample density of at least 1 point/ 
4 m 2 (AGIV, 2003) and an average density of 1 point/2 m 2 . The alti- 
metric accuracy ranges from 7 cm on a concrete surface to 20 cm 
on vegetation cover (AGIV, 2003), which is significantly better than 
the planimetric accuracy with a typical value in the order of 50 cm 
(Drosos & Farmakis, 2006). However, for less detailed and more 
common applications, the data is available in a more standard 
and practical grid format with a cell size of 5, 25 or 100 m general¬ 
ized trough Inverse Distant Weighting. The technical background 
information is summarized in Table 1. 

The original format of the LiDAR data are xyz ASCII files con¬ 
taining data points representing not only terrain elevation but 
also topographical objects like vegetation and buildings as well 
as current disturbances of the natural topography by artificial fea¬ 
tures like road banks, waste dumps, lakes and ponds. These data 
is filtered before delivery to obtain two separate datasets: one of 
the ground level points and one representing the height of topo¬ 
graphical objects like vegetation and buildings (AGIV, 2003) 
(Fig. 2). It is not the purpose of this paper to give a comprehen¬ 
sive overview of the various filter methods but summarizing, it 
can be stated that these processes are based on the significant 
differences in height values between immediate surrounding data 
points (Crutchley, 2009). This filtering however is the most criti¬ 
cal and difficult step in LiDAR data processing (Liu, 2008), as it 
creates filtering artefacts and some topographical artefacts still 
remain. Also, the additional artificial features, such as road banks 
and waste dumps, are not subject to the previously mentioned fil¬ 
ter method. For the objectives of the research project, such as 
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Fig. 1. The landscape archaeological project. (1) Study area macroscale, (2) mesoscale, (3) microscale. Known archaeological sites are indicated. 


Table 1 

Technical background information of LiDAR 
data. 

Survey data 
Co-ordinate system 

Altimetric accuracy 
Planimetric accuracy 
Minimum sample 
density 

Average sample 
density 

a AGIV (2003). 

b Drosos and Farmakis (2006). 


spatial modelling of pedogenesis and hydrology, landscape recon¬ 
struction and historical landscape analysis, a DEM is needed 
exclusively representing the natural ground level points. There¬ 
fore it is crucial to conduct a further filtering of the LiDAR ground 
points data to remove the topographical artefacts and artificial 
features. For understanding the procedure used here it is neces¬ 
sary to consider the concept of the DEM. 

In this paper the term DEM is used as described by Burrough 
(1986) in the context of “any digital representation of the contin¬ 
uous variation of the relief over space”. Three formats of digital 
models are frequently used and all possess different characteris¬ 
tics in storage efficiency and accuracy in representing the relief. 
Grid Digital Elevation Models (GDEMs) represent the relief as a 
two-dimensional regular grid whereby each grid cell contains 
one elevation value. This format is the simplest and most efficient 
in terms of storage. Consequently it allows easy manipulation but 
is relatively poor in representing the relief. Contour lines are the 
least efficient to store and manipulate data but give the best rep¬ 
resentation of the landforms. Triangulated Irregular Networks 
(TINs) represent the surface by a set of contiguous and non-over- 
lapping triangles connecting the original data points. They have 
quality properties between GDEM and contour lines (Ramirez, 
2006). 


The generation of a DEM that can be used on all scale levels is 
elaborated in five steps (Fig. 2). In the first step a GDEM is created 
based on the ground level points of the LiDAR data (dataset 2). 
This GDEM is applied to produce a filtering mask, used to filter 
the original LiDAR data from topographical artefacts and artificial 
features (step 2). Step 3 contains the filling of the empty areas 
resulting from the masking procedure. This final LiDAR dataset 
4 is used to create the final DEM representing only true ground 
points. Given the applications of this DEM on different scale lev¬ 
els, two models are derived to visualize the natural and recent 
topographical surface. First, a TIN model is created to be used 
at a large scale of analysis only (step 4). Secondly a grid model 
is produced which can be used for small- and large-scale applica¬ 
tions (step 5). 

2.2.2. Step (t) 

As the huge amount of LiDAR data is difficult to manipulate for 
filtering, an intermediate GDEM is prepared. Given that a higher 
ground resolution results in an increasing ability of recording fea¬ 
tures as well as in obtaining larger datasets, a resolution of 2 x 2 m 
is chosen in order to obtain an optimal combination of efficiency 
on the one hand and precision and accuracy of terrain representa¬ 
tion on the other hand. The high density of LiDAR data allows such 
a small grid size (Liu, 2008) that GDEM’s derived from it are most 
generally used (Ackermann, 1999; Kraus & Otepka, 2005; Liu, 
2008; Lohr, 1998; Wack & Wimmer, 2002). A TIN format is more 
used for large-scale areas because of the time consuming calcula¬ 
tions (Kraus & Otepka, 2005). The generation of a GDEM requires 
an interpolation of the irregular-spaced LiDAR points to a regular 
grid (Ackermann, 1999; Axelsson, 1999; Liu, 2008; Lloyd & Atkin¬ 
son, 2005; Wehr & Lohr, 1999a, 1999b). Despite extensive research 
and a variety of empirical work conducted by numerous authors 
(Liu, 2008), there seems to be no preferable interpolation method 
for GDEM generation (Fisher & Tate, 2006). Like several other 
authors (Axelsson, 1999; Drosos & Farmakis, 2006), a linear inter¬ 
polation method is used, based on a TIN, giving the most accurate 
and mathematically precise result. 


2001-2004 a 
Lambert BD-72/ 
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7-20 cm a 
50 cm b 
1 point/4 m 2a 

1 point/2 m 2a 
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Fig. 2. The different steps during DEM generation visualizing the natural and current topographical surface, steps (l)-(5) are explained in the text. 


2.2.2. Step (2) 

The method to extract the artificial features and topographical 
artefacts is visualized in Fig. 3. It consists of three filtering steps, 
starting with a first filtering based on the automatical creation of 
a mask using topographical vector-data. According to Liu (2008) 
“the use of additional information has the potential for increasing 
the accuracy and reliability of the filtering process”. The topo¬ 
graphical vector data chosen was surveyed at a date as close as 
possible to the acquisition date of the LiDAR data. Specific distance 
buffers are defined around these objects and used as a mask 
(Fig. 3(i)). For example, a railway has a buffer of 11 m and a small 
road a buffer of 7 m. The buffer size is based on a detailed analysis 
of the GDEM and the vector data. Given the amount and density of 
the artificial features and topographical artefacts to be filtered, it is 
recommended to use an automatical method at first. As no auto¬ 
mated filter procedure is completely accurate, additional manual 
filtering is necessary as proposed by Chen (2007) (Fig. 3(ii)). Some 
topographical artefacts and artificial features are not included in 
the automatically created mask because of an insufficient buffer 
distance or the small time gap between the LiDAR survey and the 
topographical survey on which the mask is based. To enhance 
the interpretation during this second manual filtering, the GDEM 
is visualized on screen and treated as an image. The topographical 
artefacts and artificial features are visualized maximally by histo¬ 
gram stretching and are traceable by high contrast differences in 
z-values as visible on Fig. 3(see arrow). In addition, this visualiza¬ 
tion is compared with aerial photographs and soil maps. A third 
and last filtering aims to extract small artificial elements in the 
GDEM, such as small drainage ditches and remnants of former 


ditches bordering fields. These remaining features are too small 
and too numerous to be removed manually during the second fil¬ 
tering. This third filtering implies a slope analysis followed by a re¬ 
moval of the elevation values for which the rate of the maximum 
change in z-value in a predefined radius of surrounding cells ex¬ 
ceeds a defined threshold (Fig. 3(iii)). To avoid the removal of nat¬ 
ural slopes, certain confined regions in the study area, exposing 
natural inclines exceeding this threshold, are assigned a higher 
threshold. Finally an overlay of the LiDAR dataset 2 with this final 
mask based on the three filtering steps is made. Hereby the original 
LiDAR ground points representing topographical artefacts and arti¬ 
ficial features are removed, leaving only true ground points. 

2.2.3. Step (3) 

The filtering resulted in masked areas, where new elevation val¬ 
ues are created through interpolation using two different methods. 
Small confined areas are interpolated using surrounding data 
points with a linear interpolation based on a TIN. The elevation 
for large areas is estimated based on contour lines on historical 
topographical maps of 1863 showing the landscape before the 
greater part of infrastructures was created. 

2.2.4. Step (4) 

In the fourth step, a TIN-elevation model based on the filtered 
LiDAR dataset 4 is created (Fig. 3). Unlike the grid format, the ori¬ 
ginal z-values are used which results in a more accurate modelling 
of the land surface. The TIN model fits triangular polygons to land 
facets and calculates automatically slope angle and slope orienta¬ 
tion for each triangle. To visualize the TIN, the z-values are 



Fig. 3. The different filtering steps during the removal of LiDAR points corresponding to artificial features and topographical artefacts: (i) creation of the automatic mask 
using topographical vector data, (ii) manual checking of the mask and manually adjusting when necessary, (iii) automatical filtering of ditches and field borders by slope 
analysis. 
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classified into nine classes using the equal interval method. Visual¬ 
ization is achieved by a hill shade model that enhances the micro¬ 
morphology. However, the TIN still reflects the sampling pattern of 
the LiDAR, which for example is unable to acquire data for water 
surfaces and results in a final form of artefacts. It concerns unnat¬ 
ural effects that are related to the shape and size of the triangles, 
affecting the visual impression. The filtering of incorrect data 
points from step (2) also results in similar visual effects. Therefore 
the model is interpreted as an image using ancillary data, such as 
historical map sources. Taking into account the difficulties inher¬ 
ent to the large datasets of this format, the study area is subdivided 
into a set of smaller TIN’S covering 1 km 2 resulting in an optimal 
usability for large-scale applications (Fig. 4). To meet the edge ef¬ 
fect created by this subdivision an extra border of 100 m is created, 
resulting in an overlap between adjacent TIN’S. 


2.2.5. Step (5) 

Additionally a GDEM with a cell size of 2 x 2 m is created, using 
exclusively the corrected LiDAR data points and the converted con¬ 
tour lines. Given that a TIN carries already implicitly the approxi¬ 
mate elevations of all points inside the triangles (Ramirez, 2006), 
the transformation of the randomly distributed data points to a 
GDEM is based on the previously prepared TIN. The high density 
of the LiDAR data points makes it possible to define grid cells of 
2 x 2 m. The dataset of a GDEM is much smaller than a TIN format, 
which allows small-, medium- and large-scale applications. Conse¬ 
quently, the GDEM permits a general overview of the study area. 
The dataset however, is also divided into a set of smaller grids of 
1.44 km 2 , in analogy with the TIN model (Fig. 5). This subdivision 
meets the limitations of a monitor which can only cover a range 
of 0-255 grey levels. Given the overall difference in elevation val¬ 
ues of 61 m, this results in a visualization of z-values in steps of 
0.24 m. The vertical accuracy of the source data however, lies in 
the range of centimeters resulting in an important loss of informa¬ 
tion during visualization. Subdivision into smaller areas brings 


along a significant decline of elevation range which leads to a bet¬ 
ter visualization of the topography. 

3. Results and discussion 

3.1. Interactive visualization of GDEMs 

A GDEM is not a static dataset. It can be manipulated to extract 
and visualize various types of information of the terrain surface. 
This can be done interactively. For example as represented in 
Fig. 5, the subdivision of the GDEM into 1 km 2 cells allows optimal 
visualization of the local range of z-values. Within these cells, local 
microtopography can be enhanced by a histogram stretch 
performed by spreading the presentation of the concerning height 
values over the 256 grey levels a monitor displays. This range of 
256 gray levels is however perceived by the human eye as a contin¬ 
uous gradation. When it should be possible to discriminate be¬ 
tween z-values visualized in the GDEM, not more then 10 gray 
levels should be used to represent the elevation within the display. 
Given this limitation of the human eye in distinguishing a rather re¬ 
stricted amount of colors, additional aids can help the interpreter to 
differentiate between small-scale variations of the microtopogra¬ 
phy, for example hill shading by illumination or histogram stretch¬ 
ing of the GDEM image. Especially the ability to exaggerate the z- 
values is most useful taken into account the possible enhancement 
of artefacts. To help interpret features and minimize the chances of 
features being hidden due to excessive shadow or parallel align¬ 
ment to the light source, it is advisable to create images lit at a min¬ 
imum of four different angles (Crutchley, 2009). GDEMs also offer 
the possibility to derive three dimensional images, contour lines, 
cross sections, slope models and volume calculations. 

3.2. Integrating DEM’s with historical maps and archaeological data 

Two examples of results will be discussed here. Fig. 6 shows an 
overlay of the historical map of 1909 which is the topographical 



Fig. 4. The meso-scale study area. Subdivision into (a) a mosaic of UN’s of 1 km 2 and (b) extract of the TIN elevation model derived from LiDAR data. Elevation ranges from 2.7 
to 12.8 m in this frame. Hill shading is used to enhance microtopography. Three elevation land zones can be recognized and reveal the microtopography: the low lying flat 
area of the reclaimed marshland with ditches (light gray), the sloping border (dark gray) and the convex fields on the plateau (medium gray). 
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Fig. 5. Grid Digital Elevation Model (GDEM) of the meso-scale study area with cell size of 2 x 2 m. Elevation varies from -10 to 51 m shown in 256 grayscale classes (dark is 
low; light is high). 


reference at scale 1: 20,000 just before the First World War, com¬ 
bined with an enhanced image of the TIN elevation model. The area 
shows a small elevation up to 6 m high (see contour lines), sur¬ 
rounded by roads, some even sunken roads. Clearly roads and 
houses are situated at the edge of the small bulge. The TIN clearly 
shows the microtopography with small depressions in the bulge. 
Also the field structure and field borders formed by hedge and tree 
rows follow the natural microtopography. Many fields have an arti¬ 
ficial convex topography created by raising the central parts, 
caused by digging drainage ditches around them and circular 
ploughing to the center. This field type was first reported by Snac- 
ken (1961) and archaeologically investigated by Hove (1997). Sev¬ 
eral comparable examples are showing the possibility for detailed 
interpretation of the natural conditions and the human response 
during the creation of the landscape. 

The second example illustrates how the GDEM can be used to 
detect archaeological and palaeoenvironmental features such as 



Fig. 6. Overlay of a historical map of 1909 (Military Cartographic Institute) with the 
enhanced TIN visualized as a hill shade elevation model. Elevation varies from 5.2 to 
7.0 m. 


palaeochannels or filled ditches (Fig. 7). The grid elevation model 
shows in the center an irregular depression (dark grey), some 
0.3 m lower than the surrounding grounds (light grey-white). A 
channel structure is marked in the middle. This pattern represents 
a north-south orientated palaeochannel or filled ditch which 
clearly can be traced in the GDEM over a distance of 2.5 km. Other 
data sources (Fig. 7: nrs. 2-7) give no evidence at all of this feature. 
These findings made a more purposive localization of auger tran- 
sections for geophysical, palaeoecological, geomorphological and 
archaeological survey possible. 

Assembling the GDEM and the TIN with historical maps and 
archaeological data forms the basis of further spatio-temporal 
landscape change analysis of the study area. 

4. Conclusions 

The aim of this paper was to generate a high precision DEM, 
representing the natural and recent topography and to review 
the possibilities of this DEM for historical geographical purpose. 
Trough the collection of elevation points at a much higher density 
than traditional survey methods, high density airborne LiDAR data 
forms the basis for this high quality DEM, meanwhile reducing 
data volume and increasing efficiency. As only true ground points 
should be implemented, three filtering steps were performed start¬ 
ing with an automatic removal based on topographical vector data 
followed by further manual processing by visual interpretations 
and an automatic filtering using slope analysis. Where artificial 
features and artefacts were extracted, ground heights were pre¬ 
dicted based on a linear interpolation method using surrounding 
data points and historical topographical maps of 1863. Based on 
this set of filtered LiDAR data two different models were created. 
At first, a TIN is derived, which is used for large-scale applications 
and results in the best fitting to the terrain. The second model cre¬ 
ated is the GDEM with a pixel size of 2 by 2 m, being most efficient 
in terms of storage and manipulation. This grid will be used for 
small-, medium- and large-scale applications. Both DEM models 
represent the natural topography free of artefacts caused by mod¬ 
ern objects and infrastructures created since the 19th century. De¬ 
spite the considerable expenses for data acquisition and enormous 
volume of data, LiDAR is highly recommended for terrain visualiza¬ 
tion trough DEM generation. LiDAR derived DEMs offer the oppor¬ 
tunity to study the earth’s topography in such a high resolution 
which no other survey and data source can deliver for similar 
small-scale research. Although it is a semi-continuous model of a 
continuously variable being the terrain surface, it is a valuable 
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Fig. 7. Composite of the GDEM with thematic and historical maps. (1) GDEM 2004, (2) cadastral map of 1576, (3) soil map of 1960, (4) geomorphological map of 1999, (5) 
historical map of 1576, (6) historical map of 1775 and (7) historical map of 1910. 


complementary source for historical landscape research providing 
valuable, unique and high resolution information. The analysis of 
LiDAR derived DEMs requires a similar approach as aerial photo¬ 
graphs interpretation. The interpretation of LiDAR based DEMs 
depends highly on the visualization technique and the resolution 
and should be manipulated interactive in order to detect small 
topographical details, for example trough hill shading. Within the 
scope of this interdisciplinary project the results of the DEM gener¬ 
ation will be integrated with the outcome of the research of the 
natural sciences and the archaeological team with the aim of relat¬ 
ing the settlement history and the main land use changes to natu¬ 
ral environmental conditions and the changes of landscape 
ecological processes. 

There can be concluded that the optimal use of LiDAR data de¬ 
mands multi-scale processing and the use of both grid and TIN- 
formats. LiDAR derived DEMs allow a more detailed landscape 
analysis of microtopography, field-, settlement- and road patterns 
and become a common reference data layer for various 
disciplines. 

Acknowledgements 

The realization of this project, especially the acquisition of the 
expensive LiDAR data would not be possible without the funding 
outgoing from the University of Ghent. More information about 
the project can be found on http://cartogis.ugent.be/sandy_flan- 
ders/. Also our special thanks to the Agency for Geographic Infor¬ 
mation in Flanders (AGIV) and the National Geographical 
Institute (NGI-B) whose willing cooperation in assembling various 
source materials was essential. The authors also want to thank 
Msc. Iason Jongepier for the acquisition of historical documents 


and Dr. Prof. A. De Wulf for the technical support concerning the 
processing of LiDAR data. 

References 

Ackermann, F. (1999). Airborne laser scanning - present status and future 
expectations. ISPRS Journal of Photogrammetry and Remote Sensing, 54, 64-67. 

AGIV. (2003). Digitaal Hoogtemodel Vlaanderen. Nieuwsbrief GIS Vlaanderen 
(pp. 3-21). 

Antrop, M., Eetvelde, V. V. (2007). L’atlas des paysages flamands dans la politique 
integree d’amenagement dut territoire. In Berlan-Darque, M., Luginbiihl, Y., 
Terrasson, D. (Eds.), Paysages: de la connaissance a I’action. Editions Quae, c/o Inra, 
Collection Update Sciences E? Technologies (pp. 139-148). Versailles. 

Axelsson, P. (1999). Processing of laser scanner data - algorithms and applications. 
ISPRS Journal of Photogrammetry and Remote Sensing, 54, 138-147. 

Bender, O. (2009). The concept of a historic landscape analysis using GIS with focus 
on Central Europe. In Bender, O., Evelpidou, N., Krek, A., Vassilopoulos, A. (Eds.), 
Geoinformation technologies for geocultural landscapes, European perspectives (pp. 
129-144). London. 

Bewley, R. H., Crutchley, S. R., & Shell, C. A. (2005). New light on an ancient 
landscape: lidar survey in the Stonehenge World Heritage Site. Antiquity, 79, 
636-647. 

Burrough, P. A. (1986). Principles of geographical information systems for land 
resources assessment. New York: Oxford University Press. 

Carey, C., Brown, T., Challis, K., Howard, A., & Cooper, L. (2006). Predictive modelling 
of multiperiod geoarchaeological resources at a river confluence: A case study 
from the Trent-Soar, UK. Archaeological Prospecdon, 13, 241-250. 

Challis, K. (2006). Airborne laser altimetry in alluviated landscapes. Archaeological 
Prospecdon, 13, 103-127. 

Challis, K., Kokalj, Z., Kincey, M., Moscrop, D., & Howard, A. J. (2008). Airhorne lidar 
and historic environment. Antiquity, 82, 1055-1064. 

Chen, Q.. (2007). Airborne LiDAR data processing and information extraction. 
Photogrammetric Engineering and Remote Sensing, 73, 109-112. 

Crutchley, S. (2009). Shedding light on the past: Using airborne LIDAR to understand 
ancient landscapes. In Bender, O., Evelpidou, N., Krek, A., Vassilopoulos, A. (Eds.), 
Geoinformadon technologies for geocultural landscapes, European perspectives (pp. 
87-106). London. 

Devereux, B., Amable, G., & Crow, P. (2008). Visualization of LiDAR terrain models 
for archaeological feature detection. Antiquity, 82, 470-479. 
























8185 


/. Werbrouck et al.f Expert Systems with Applications 38 (201 2 ) 8178-8185 


Drosos, V., Farmakis, D. (2006). Airborne laser scanning and DTM generation. 
In Manolas, E. (Eds.), Naxos international conference on sustainable 
management and development of mountainous and island areas (pp. 206- 
218). Greece. 

Eeckhaut, M. V. D., Verstraeten, G., 82 Poesen, J. (2007). Morphology and internal 
structure of a dormant landslide in a hilly area: The Collinabos landslide 
(Belgium). Geomorphology, 89, 258-273. 

English Heritage Aerial Survey and Investigation English Heritage. Available from 
http://www.english-heritage.org.uk/server/show/nav.8730 (1.09.2009). 

Fernandes, A. M., Utkin, A. B., Lavrov, A. V., 82 Vilar, R. M. (2004). Development of 
neural network committee machines for automatic forest detection using lidar. 
Pattern Recognition, 37, 2039-2047. 

Fisher, P. F., & Tate, N. J. (2006). Causes and consequences of error in digital 
elevation models. Progress in Physical Geography, 30, 467-489. 

Gallagher, J., & Josephs, R. (2008). Using LiDAR to detect cultural resources in a 
forested environment: An Example from Isle Royale National Park, Michigan, 
USA . Archaeological Prospection, 15, 187-206. 

Harmon, J. M., Leone, M. P., Prince, S. D., & Snyder, M. (2006). Lidar for 
archaeological landscape analysis: A case study of two eighteenth-century 
Maryland plantation sites. American Antiquity, 71, 649-670. 

Heyse, I. (1983). Preliminary results of the study of a vistulian late glacial drainage 
pattern in the scheldtbassin (Belgium - Flemish Valley - Moervaart depression). 
Quaternary Studies in Poland, 4, 135-143. 

Hove, R. V. (1997). De “klassieke” bolle akkers van het Waasland in archeologisch 
perspectief Sint-Niklaas. 

Kasai, M., Ikeda, M., Asahina, T., & Fujisawa, K. (2009). LiDAR-derived DEM 
evaluation of deep-seated landslides in a steep and rocky region of Japan, 
LiDAR-derived DEM evaluation of deep-seated landslides in a steep and rocky 
region of Japan. Geomorphology, 113, 57-69. 

Kraus, K., & Otepka, J. (2005). DTM modelling and visualization - the SCOP 
approach. Photogrammetric Week 05 Wichmann Verlag, 241-252. 

Liu, X. (2008). Airborne LiDAR for DEM generation: Some critical issues. Progress in 
Physical Geography, 32, 31-49. 


Lloyd, C. D., & Atkinson, P. M. (2005). Deriving ground surface digital elevation 
models from LiDAR data with geostatistics. International Journal of Geographical 
Information Science, 20, 535-563. 

Lohr, U. (1998). Digital elevation models by laser scanning. Photogrammetric Record, 
16, 105-109. 

Mutlu, M., Popescu, S. C., 82 Zhao, K. (2008). Sensitivity analysis of fire behavior 
modeling with LiDAR-derived surface fuel maps. Forest Ecology and 
Management, 256, 289-294. 

Powlesland, D., Lyall, J., Hopkinson, G., Donoghue, D., Beck, M., Harte, A., et al. 
(2006). Beneath the sand - remote sensing, archaeology, aggregates and 
sustainability: A case study from Heslerton, the Vale of Pickering, North 
Yorkshire, UK. Archaeological Prospection, 13, 291-299. 

Ramirez, J. R. (2006). A new approach to relief representation. Surveying and Land 
Information Science, 66, 19-25. 

Saye, S. E., Wal, D. v. d., Pye, K., 82 Blott, S. J. (2005). Beach-dune morphological 
relationships and erosion/accretion: An investigation at five sites in England 
and Wales using LIDAR data. Geomorphology, 72, 128-155. 

Snacken, F. (1961). Streekindeling en begrenzing van het Land van Waas. Belgisch 
Tijdschrift van de Vereniging voor Aardrijkskundige Studies, 30, 217-255. 

Utkin, A. B., Lavrov, A. V., Costa, L., Simoes, F., 82 Vilar, R. (2002). Detection of small 
forest fires by LIDAR. Applied Physics, 74, 77-83. 

Wack, R., 82 Wimmer, A. (2002). Digital terrain models from airborne laser scanner 
data - a grid based approach. International Archives of Photogrammetry, Remote 
Sensing and Spatial Information Sciences, 34, 293-296. 

Wehr, U., 82 Lohr, U. (1999a). Airborne laser scanning - an introduction and 
overview. ISPRS Journal of Photogrammetry and Remote Sensing, 54, 68-82. 

Wehr, U., & Lohr, U. (1999b). Theme issue on airborne laser scanning. ISPRS Journal 
of Photogrammetry and Remote Sensing, 54, 61-63. 

Zhou, G., Song, C., Simmers, J., & Cheng, P. (2004). Urban 3D GIS from LiDAR and 
digital aerial images. Computers Ef Geosciences, 30, 345-353. 

Zhou, G., 82 Xie, M. (2009). Coastal 3-D morphological change analysis using LiDAR 
series data: A case study of Assateague Island National Seashore. Journal of 
Coastal Research, 25, 435-447. 


